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n-lmpulse, Gravity-Assist Trajectory Optimization Q ***** 


Optimal number of maneuvers not 
known a priori 

- Deep-space maneuvers (DSMs) frequently 
improve performance 

- Two (or more) DSMs can be optimal for 
tightly/uniquely constrained trajectory legs 


Optimal EEJ trajectory w / 2 DSMs on 1 leg 


lee 


Event # 6: 
insertion 
Jupiter 
3/3/2027 

— 5. 59 3 km/ » 
DEC = -7.0 
Ar = 1,284 feiri/s 

rn = 670 kg 


Missions interested in best performance 
possible, i.e., the global optimum 

- Can enable cost- or time-constrained missions 

- Often interested in more than 1 objective, 
e.g., max. delivered mass & min. TOF 

Problem challenges: 


Event # 3i 
chemical burn 
deep-space 
7/12/2024 
Ar = 2,230 km/ 3 
m = 1S35 kg 


x (km) 
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Event # 2; 

chemical burn 

deep-space 

1/5/2024 

Air = 0.950 km/ s 

m = 3734 kg 


Event # 5; 
chemical bum 
deep-space 
9/15/2024 
Ait = 1.050 km/* 
m — 1309 kg 



Event # 4: 
un powered flyby 
Earth 
9/6/2024 

= 9,080 km/ s 
DEC - 1.2^ 
altitude = 300 km 
m — 1835 kg 


E 
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Event# l: 
launch 
Earth 
9/3/2023 

C3 = 10.066 ft m s /* 2 
DLA = 8.4 
m = 5054 kg 
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- Grid searches can be intractable for any number & variety of DSMs 

- Optimization requires initial guess 

- Maneuvers and gravity assists can create highly sensitive optimization problems 
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Objective 


a 




Globally optimize chemical propulsion based trajectories 
with an arbitrary number of maneuvers & gravity assists 


Method should be: 

• Automated 

• No requirement of a user-defined initial guess 

• Able to search broad design space 

• Efficient (medium-fidelity is appropriate) 

• Capable of handling multiple objectives 
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Prior Approaches to Global Search 


Q atatSA 


• Grid search 

- Lambert scans over range of event dates & flyby 
bodies 

- Strategies developed to include one maneuver per leg 

• Patel and Longuski -> STOUR (Purdue, JPL) 

• Lantukh and Russell (UT-Austin) 

• Maneuvers limited to specific type v-infinity leveraging or 
broken plane maneuver 


• Stochastic search 

- Strategically sample design space 

- Typically use a direct, Lambert-based trajectory formulation: 

• Vary departure date, time to DSM, & DSM components 

• Propagate forward to DSM point 

• Solve Lambert problem to subsequent body, repeat for all legs 

• Guaranteed feasibility for unconstrained problems 

- Not limited in maneuver type, but problems are frequently very sensitive to initial guess 

- Most approaches not capable of more than 1 DSM 

- Global-local hybrid scheme can improve efficacy 


Example STOUR plot (TOF vs. LD) 

Petropoulos et al., “Trajectories to Jupiter...” JSR, 2000 

PATH: 3 2 3 3 5 
VINF(KM/S): 3.00 3.25 3.50 3.75 



LAUNCH DATES SEARCHED: 99/ 1/ 1 TO 32/ 1/ 1 by 5.0 days 
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MG AnDSMs Trajectory Transcription 
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Multiple gravity assist with n deep space maneuvers using shooting scheme s 
MGAnDSMs 

- Aim for robustness & efficiency with an arbitrary number of DSMs with a direct formulation 

- Nonlinear programming problem 

- Employs forward/backward shooting similar to Byrnes & Bright (CATO) 

- Nominally Kepler propagation between maneuvers 

- Analytic match point (MP) constraints 

- User specifies n, optimizer reduces 1 or more DSM magnitude to zero if <n optimal 


n = 4 maneuvers 


Right MP state: 



Phase 

End 

Point 


Core Optimization 

Variables 


At, 


phase; 


Ak 


At,, ...At, 


At, 


n+1 


Av 1 ; Av 2 , Av n 

^initial 

^final 


Match point 

continuity constraint. c mp \?x k > fy ty > k k > kc > ^ y ^ z z > rrr ] £ 


5 




MGAnDSMs Mission Formulation 
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• Multiple gravity assists accommodated in mission structure 

• Allows for variation of flyby bodies at control nodes of transcription 

- Journeys start & end a user-required bodies/states 

- Journeys can be composed of multiple phases with variable flyby bodies to improve performance (“null 
gene” approach) 

• Zero sphere of influence patched conics 

- Flyby constraints ensure physically realizable gravity assist 

• Flexible to numerous mission constraints 
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MGAnDSMs Global Optimization 
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Combine monotonic basin hopping (MBH) & sequential quadratic programming 
(SQP) MBH+SQP 

Stochastic, global search scheme 
No initial guess required 

Adept at multi-modal problems w/ clustered local minima 
Stochastic “hops” evaluated from base solution 




Multi-objective Optimization 
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ts 


• MGAnDSMs structured within multi-objective hybrid optimal control 
algorithm (discrete & continuous variables) 

• Multi-objective genetic algorithm (GA) serves as outer loop systems 
optimizer around direct-method inner loop trajectory optimizer 

- Outer loop: non-dominated Sorting Genetic Algorithm II (NSGA-II) searches over discrete 
mission parameters, defining trajectory problem for inner loop 

• Variables include: Flyby body, target body, launch vehicle, launch C3, launch epoch 

- Inner loop: MBH+SQP solves trajectory problem & establishes obj. func. values 

- Generates representation of Pareto front (optimal tradeoff between objectives) 


Initial generation 


Final generation 
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Delivered 

Mass 



*■ 


Population evolves via 
genetic operators Delivered 

^ Mass 
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TOF 


TOF 
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Cassini Test Case Comparison 
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• Established global optimum serves a test case for comparison of MGAnDSMs & a Lambert- 
based transcription (MGADSM/c); one DSM allowed 

• MGAnDSMs identified global optimum in 10 out of 10 cases 

• MGADSk identified solution w/in 50 m/s of global optimum, but not the global optimum (1 .004 
km/s) after an 8 hour run time 

• Order of magnitude improvement in time to best solution 
- Median time is 25x faster 


Run 

Number 

MGAn DSM s 
Time to Identify 
Global Optimum 
(minutes) 

MGADSM/r 
Time to Identify 
Best Solution 
(minutes) 

MGADSM/r 

Time Identify Solution 
within 50 m/s of 
Global Optimum 
(minutes) 

MGADSM/c 
Best AV 
(km/s) 

1 

5.3 

223.8 

220.3 

1.025 

2 

7.9 

140.9 

32.9 

1.043 

3 

8.8 

260.8 

149.3 

1.026 

4 

12.8 

74.3 

74.3 

1.014 

5 

63.1 

178.1 

122.0 

1.020 

6 

4.3 

295.4 

112.0 

1.034 

7 

5.0 

471.3 

465.7 

1.039 

8 

6.3 

382.9 

328.5 

1.049 

9 

9.4 

80.6 

58.4 

1.033 

10 

38.4 

63.6 

63.6 

1.038 

Mean 

16.1 

217.2 

162.7 

1.032 

Median 

8.4 

200.9 

117.0 

1.033 


VVEJ 



9 




Example Problem: Multi-objective Optimization to 
Jupiter 
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• Evaluated MGAnDSMs multi-objective capability on a multiple gravity assist to Jupiter with up to 
2 DSMs per phase 

• Flyby bodies are varied with any combination up to 5 bodies 

• Three objectives: maximize log 10 (final mass), minimize TOF, minimize Jupiter arrival C3 

- Pareto front: 3D surface of equally optimal solutions 

• Population size of 256 and the MBH+SQP inner loop is allowed to run for 40 minutes 


Outer-loop Design Variables Common Mission Parameters 


Design Variable 

Value 

Resolution 

Launch window 
open epoch 

{1/1/2021,1/1/2022, 
1/1/2023, 1/1/2024} 

1 year 

Flyby body 

{Venus, Earth, Mars, 
null, null, null} 

n/a 

Flight time 

[730, 3467.5] days 

182.5 days 


Description 

Value 

Launch window 

365.24 days 

Launch declination 

[-28.5, 28.5] deg 

Launch vehicle curve 

Atlas V, 551 

Chemical l sp 

320 s 

Jupiter arrival date 

Determined by 
optimizer 

Jupiter insertion orbit semi-major axis 

10,054,900 km 
(140.6 Rj) 

Jupiter insertion orbit eccentricity 

0.911 

Maximum number of DSMs 

2 

Inner-loop objective function 

Max: log 10 (final mass) 

Inner-loop run time 

40 minutes 
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Delivered mass (kg) 


Best Non-Dominated Front 


• Representation of 3D Pareto front generated after 1 00 generations 

• 3.5 days of run time on 64-core processor 


2500 

2000 


1000 


0 
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Arrival v-infinity (km/s) 



Example Optimal Jupiter 


Highest Delivered Mass Trajectory 

• EVEEJ sequence 

• Delivered mass: 31 92 kg 

• TOF: 7.9 years 

• Arrival C3: 34.1 km/s 



y (km) 


jectories 
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Shortest TOF Trajectory 

• Direct EJ (2 DSMs optimal) 

• Delivered mass: 642 kg 

• TOF: 2.0 years 

• Arrival C3: 45.1 km/s 
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y (km) 



Optimal Jupiter Trajectory with 2 DSMs in 1 Phase Q 


EVEJ sequence 
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Conclusions 


Q 


• Developed high-thrust trajectory optimization transcription for any number of 
maneuvers between flybys 

• Multiple shooting framework & analytic derivatives provide robustness and 
enable outer-loop efficiency 

• Formulation allows for an automated, global search without a user-supplied 
initial guess 

• Capability to generate Pareto-optimal solutions using multi-objective hybrid 
optimal control algorithm 

• General applicability to almost any variety of interplanetary mission 

- Flexible to unique trajectory/mission constraints 

• Large problems can become computationally tractable 
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Backup 
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Multi-objective Optimization 
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• Want to optimize any number of mission design metrics 

- e.g., payload mass, TOF, arrival C3 

- Often coupled & competing 

- Fully map mission trade-offs between optimal solutions 


• Optimize multiple objectives 
simultaneously 

- Entire set of optimal solutions 

- Goal: generate representation of 
Pareto front 

- Traditionally use repetitions of 
single objective technique 
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NSGA-II 
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• Develops globally-optimal Pareto solutions using non-dominated sorting 

- Conducts stochastic, global search with population of designs 

• Fitness assignment based on “nearness” to Pareto front 


- x 1 dominates x 0 if: 

Vp:/ P (x,)</ P (x 2 ) 

p = l,2,...fl obj 


and 

3 p:/ P (Xi)</ p (x 2 ) 

p 1,2 

C/) 

c n 

If neither design dominates other, they 

Cu 

E 
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are non-dominant 


H — 


• Non-dominated sorting: 


- Assign fitness based on design’s non- 
dominated front 

- Designs closer to Pareto front get more 
mating opportunities 
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Genetic Algorithm 


Initial random population 
generated 
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Models Darwinian evolution 
- Mimic natural selection & reproduction 
Searches with population of designs 
Globally search design space 
No initial guess required 


18 



